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ABSTRACT 



Q^ I The question of whether flares occur as a Poisson process has important consequences 

for flare physics. Recently Lepreti et al. presented evidence for local departure from 
^ ■ Poisson statistics in the Geostationary Operational Environmental Satellite (GOES) X- 

^^ ■ ray flare catalog. Here it is argued that this effect arises from a selection effect inherent 

in the soft X-ray observations; namely that the slow decay of enhanced flux following 
a large flare makes detection of subsequent flares less likely. It is also shown that the 
power-law tail of the GOES waiting-time distribution varies with the solar cycle. This 

^^ \ counts against any intrinsic significance to the appearance of a power law, or to the 

i^ ' value of its index. 

^ ■ Subject headings: Sun: activity — Sun: flares - Sun: corona - Sun: X-rays 

-t— > 



Introduction 



X 

j^ ■ There has been increased recognition of the importance of flare statistics for understanding the 

mechanisms of energy storage and release in the solar corona. The distribution of times At between 
events ("waiting times") provides information on whether flares occur as independent events, and 
also represents a test for certain flare models. In particular, the avalanche model for flares (Lu and 
Hamilton 1991; Lu et al. 1993) and variants thereof, predict that flares are independent and occur 
with a constant rate, and hence obey Poisson statistics. The waiting-time distribution (WTD) for 
a Poisson process with a rate A is an exponential, 

P(At) = Ae'^^*. (1) 

Observational determinations of the flare WTD based on different data sets have given a variety 
of results (see e.g. Wheatland et al. 1998 for a brief review). Boffeta et al. (1999) showed that the 
WTD constructed from the GOES catalog of 20 years of flaring exhibits a power-law tail. They 
argued that the appearance of a power law contradicts the avalanche model prediction of Poisson 



statistics, and supports instead an MHD turbulence model of the flaring process. Special significance 
was attributed to the value of the power-law index (a ~ —2.4) derived from the observations. In 
response, Wheatland (2000) argued that over the course of the several solar cycles included in the 
GOES data, the flaring rate varies by more than an order of magnitude, and so flares cannot be 
assumed to be occurring at a constant rate. If the flaring process can be represented by a piecewise 
constant Poissson process with constant rates Aj for intervals tj, then the WTD may be described 
by 

P{At) = Y,¥^i^ie-^'^\ (2) 

i 

where cpi = Xiti/ J2i ^i^i is the fraction of events corresponding to the rate Aj. Rates Aj and intervals 
ti were estimated for the GOES catalog events using a Bayesian procedure. The resulting WTD 
[equation (2)] was shown to qualitatively reproduce the observed power-law tail. The GOES catalog 
involves flares occurring in all active regions on the Sun. Wheatland (2001) extended the analysis 
to flares in individual active regions, and showed that equation (2) accounts quantitatively for the 
observed WTDs in a number of very flare-productive active regions. 

Recently Lepreti et al. (2001) responded by showing that a statistical test rejects the hypothesis 
that the GOES catalog events are locally Poisson. They showed that the flare WTD can be fitted 
by a Levy function, and argued that this reflects the existence of memory in the system and once 
again supports a turbulence model for the flare process. 

Although the question of whether flares occur as a Poisson process may seem arcane, the 
consequences for the understanding of the flare process are substantial. 

In this paper the result of Lepreti et al. (2001) is re-examined. In § 2 an alternative explanation 
for the local departure of the GOES flares from Poisson statistics is given, namely that it arises 
from the failure to detect flares occurring soon after large flares because of the increase in soft 
X-ray flux associated with the large flare. The nature of soft X-ray data and the event selection 
procedure used to compile the GOES catalog makes this effect inevitable. In § 3 it is shown that 
the power-law tail of the GOES WTD varies with the solar cycle. This result argues against any 
particular significance to the value of the power-law index of the tail of the WTD. Finally in § 4 
the consequences of these results to our understanding of flares are discussed. 



2. Local departure from Poisson statistics 

Lepreti et al. (2001) applied a test for local Poisson behaviour devised by Bi et al. (1989) to the 
GOES catalog events. Flares are replaced by their peak times. Let Xi denote the time interval 
between the peak time of event i and that of its nearest neighbour. Also, let Yi denote the time 
interval to the other neighbour. Under the local Poisson hypothesis, Xi and Yi have probability 



distributions P{Xi) = 2Ai exp(— 2AjXj) and P{Yi) = Aj exp(— Ajl^). In this case the variable 

Xi + 2^ 

should be uniformly distributed between zero and one, and then the observed distribution of Hi 
provides a test of the local Poisson hypothesis. Lepreti et al. (2001) reported a significantly non- 
uniform distribution of Hi for the GOES events. We have repeated the procedure for GOES events 
between 1981-1999 above CI class (for details of the GOES data, see Wheatland 2001), and the 
result is shown in Fig. 1. The solid curve is the cumulative distribution oi Hi. A Kolmogorov- 
Smirnov test rejects the hypothesis that the distribution of Hi is uniform, confirming the finding 
of Lepreti et al. (2001). 

Of course, the departure of the GOES catalog events from Poisson behaviour does not neces- 
sarily mean that solar flares are not locally Poisson. It is important to consider how the catalog was 
constructed, and in particular any selection biases that might introduce a departure from Poisson 
statistics. Recently Wheatland (2001) identified such an effect, which was termed obscuration. The 
GOES detectors record soft X-ray light curves for flares, which are characterised by relatively fast 
rises to a peak flux, followed by slow decays in flux due to the cooling of hot plasma in the corona 
produced by the flare (Feldman 1996). Large flares produce enhancements of more than a factor 
of 100 in flux, and this increase may take many hours to decline. During the decay time subsequent 
flares may be missed because of the enhanced background. The selection procedure used to compile 
the GOES catalog imposes constraints on the detection of flares. The start of a flare is defined by 
four one-minute flux values that are monotonically increasing, with the final flux being at least 1.4 
times the first. For flares occurring shortly after a large flare, this implies that the second flare must 
produce an increase in flux at least of order 40% of the flux of the first flare to be detected. Even 
quite significant events can be missed in the wake of a large flare. Since flares follow a power-law 
peak flux distribution (Hudson 1991), the majority of flares are small and a large number of flares 
are expected to be missing from the catalog as a result. Wheatland (2001) estimate that in the 
absence of obscuration the number of flares above CI class would be higher by (75 it 23)%. 

The following simulation conflrms that obscuration can account for the observed departure 
from Poisson statistics. Using the Bayesian decomposition of the GOES catalog into rates Aj and 
intervals tj described in Wheatland (2000), a sequence of flare times was generated as a piecewise 
constant Poisson process. As a check, the Bi et al. (1989) test was applied to this sequence of 
times, and local Poisson behaviour was conflrmed (the cumulative distribution of Hi is shown by 
the dotted line in Fig. 1). Next, a decay time was assigned to each synthetic flare, based on 
random selection from the decay times observed for the real GOES events. To mimic obscuration, 
any events falling within the decay time of another event were excluded from the sequence. The Bi 
et al. (1989) test was applied to the resulting sequence of flares. The result is shown as the dashed 
curve in Fig. 1. Once again the distribution of Hi is signiflcantly different from that expected from 
local Poisson behaviour. The distribution is clearly very similar to that found for the GOES data. 
Although this simulation is crude, it shows that a Poisson simulation with obscuration produces 



an effect mimicking that reported by Lepreti et al. (2001) for the real GOES flares. 



3. Solar cycle variation 

Boffeta et al. (1999) argued that the observed power-law tail in the WTD for the GOES events 
is "robust," and even claimed that the power law in waiting times is "as firmly established as the 
power laws observed for total energy, peak luminosity, and time duration." Lepreti et al. (2001) 
went further by assigning special significance to the fitting of the power law by a Levy function. 
The value of the power-law index is important to this fit. 

A feature of the other power laws (energy, peak fiux, and possibly duration) observed for solar 
flares is that they do not vary substantially with the 11-year solar cycle (Dennis 1985; Lu and 
Hamilton 1991), although evidence has been presented for small variations (Bai 1993). Here we 
show that, by contrast, the power law in waiting times is strongly dependent on the phase of the 
cycle. 

Fig. 2 shows the WTDs for the maximum and minimum phases of the solar cycle. Maximum 
phase has been defined as times when the smoothed monthly sunspot number^ is above 125, and 
minimum has been defined as times when the sunspot number is below 30. These definitions are 
arbitrary, but they provide an objective choice of cycle phase, i.e. they are made without reference 
to the fiare data. The two histograms are clearly different. The histogram extending to larger 
waiting times corresponds to minimum phase. This is not surprising: the mean rate of flaring 
varies by more than an order of magnitude over the cycle, and since the mean rate is the reciprocal 
of the average waiting time, it follows that the WTD must vary with the cycle. Fits to the power- 
law tails (> 10 h) of the WTDs shown in Fig. 2 give power-law indices of —3.2 it 0.2 and — 1.4zb 0.1 
for the maximum and minimum phases respectively. The solid curves shown in Fig. 2 correspond to 
the piecewise constant Poisson model (2) with rates and times obtained by the Bayesian procedure 
(Wheatland 2000). This model accounts for the qualitative shape of the WTDs in each case, 
although clearly there is some discrepancy, which is most likely due to the failure of the Bayesian 
procedure to account for all of the rate variation. 



4. Conclusions 

Lepreti et al. (2001) demonstrated a departure from Poisson statistics in the GOES catalog of flares, 
and attributed this to "time invariance" in the flaring history. Here an alternative explanation is 
given, namely that the flux increase due to large flares reduces the likelihood of detecting subsequent 
smaller flares. This selection effect (termed obscuration) introduces a departure from Poisson 



^This data, as well as the GOES catalog, is available from ftp://ftp.ngdc.noaa.gov/STP/SOLAR_DATA. 



statistics because it effectively means that the rate decreases when a large flare occurs, and hence 
the flares in the GOES catalog are not strictly independent. A simulation conflrms that obscuration 
produces an effect very similar to that reported by Lepreti et al. (2001). This study underlines the 
need to consider all potential biases in the statistical study of flares. 

The power-law tail to the WTD constructed from the GOES catalog has been shown to vary 
with time, and in particular with the solar cycle. On this basis it is difficult to argue that the power 
law is "robust," in marked contrast with the other power laws reported in flare statistics. 

The picture emerging from this study and others (e.g. Wheatland 2001) is that flares occur as 
a time-varying, or non-stationary Poisson process. The observed WTD depends on the observed 
rates of flaring (Wheatland 2000) which naturally accounts for the variation of the WTD over 
the solar cycle. There is no particular significance to the appearance of a power law in the WTD. 
Indeed, adopting the continuous version of (2), viz. 

1 1"°° 
P{At) = - /(A)AV^^*ciA, (4) 

-^0 Jo 

where f{X)dX is the fraction of time that the rate is in the range A to A -|- dA and Aq = /g A/(A) dX 
is the mean rate of fiaring, the asymptotic behaviour of P{At) is obtained by Taylor expanding 
/(A) about A = 0. This gives 

P(At) = M2) (At)-3 + ^4M(At)-4 + ... (5) 

Ao Ao 

Equation (5) suggests that the asymptotic form of P{At) will always be a power law, with an index 
of —3. In assessing this result it should be remembered that (4) is an approximate expression (it 
assumes the rate does not varying greatly during a waiting time). Lepreti et al. (2001) appeared to 
forget this point in criticising the WTD obtained by Wheatland (2000) with an exponential form 
in (4). 

Interestingly, Earthquake data also often exhibit power-law WTDs due to non-stationarity of 
Earthquake sequences. Utsu (1970) states that there is "no primary importance" to the appearance 
of a power law, but the power law exhibited by Earthquake size distributions (the Gutenberg- Richter 
relationship) is considered important. The situation seems to be closely analogous to that for fiares. 

The results of this paper appear to be consistent with an avalanche model of fiaring subject to 
a time- varying rate of driving (Norman et al. 2001). However, the appearance of Poisson statistics 
does not exclude other models with similar general hypotheses concerning energy storage. Indeed, 
any steady state system with rates of downward transition in energy (fiares) depending only on 
the size of the transition, and with a mean energy larger than the majority of flares, will exhibit 
Poisson statistics (Craig and Wheatland 2001). 
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Fig. 1. — The Bi et al. (1989) test for Poisson behaviour, applied to the GOES data (sohd), to a 
Poisson simulation (dotted), and to the same Poisson simulation including obscuration (dashed). 
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Fig. 2. — WTDs for the maximum and minimum phases of the solar cycle. 



